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Abstract. We discuss the consequences of the spectral and morphological modification of galaxy cluster radio 
halos due to the Sunyaev-Zeldovich (SZ) effect for the interpretation of existing and upcoming high frequency 
radio observations. Likely these modifications have affected the interpretation of the existing Goma cluster radio 
data. The radio halo emission visible at low 5 GHz) frequencies is at higher (> 10 GHz) frequencies completely 
over-compensated by the thermal SZ decrement. Thus, the total radio emission of a galaxy cluster goes through 
zero (in comparison to the constant cosmic microwave background (CMB) emission) at a frequency of several 
GHz. Since the radio halo brightness has a narrow radial profile compared to the SZ decrement, a central emission 
region is surrounded by a decrement within the intermediate frequency range of several GHz. The size of this 
emission regions shrinks with increasing frequency until the decrement dominates everywhere in the cluster. 

Key words. Intergalactic medium - Galaxies: cluster: general - Radio continuum: general - Radiation mechanism: 
thermal - Radiation mechanism: non-thermal - Gosmic microwave background 


1. Introduction 


also the self-annihilation of neutralino dark matter par¬ 


ticles was pro posed as an electron source ( |Colafrancesco 


plasma. Some - if not all - clusters host ultra-relativistic 
particle populations, observable by synchrotron emission 
of ~ 10 GeV electrons. This emission forms the cluster 
radio halos. Recently, the number of detected cluster ra- 


Cluste. s of galaxies are hlled by a hot (fcJe ^ b keV) ^ Mele| 12001 ) although the observed correlation of radio 


dio halos increased significantly (Giovannini et al. 

1999; 

iGiovannini & Feretti 

O 

o 

o 

Kempner & Sarazin 

pool 

) and 


the prospects for upcoming sensitive radio telescopes to 
detect large numbers of cluster radio halos are promising 


(EnBlii |L fc RoLLgerhi g-2Q02) 


halos with merging events clearly disfavours this scenario. 
Nevertheless, such dark matter based electron sources may 
lead to a lower level of radio emission and a search for 
such emission has the potential to constrain scenarios of 
the nature of the dark matter ( Bohm et al. 2002 ). 

Such desirable high frequency measurements of the 
cluster halo emission will unavoidably be affected by the 
thermal SZ effect ([Zeldovich fc Sunyaev 1969; Sunyaev & 


Radio halos seem to be triggered by strong cluster 


merger ( Buote 2001 ), although the detailed physical mech¬ 
anism providing the energetic electrons is still to be re¬ 
vealed. Two classes of theoretical models are discussed in 
the literature: i) in-situ acceleration of the rapidly cool¬ 
ing electrons by tur bulence an d shock waves throughout 
the cluster volume ( Jaffe 1977 , and others), and ii) sec¬ 
ondary electron injection after particle physical interac¬ 
tions. In-situ acceleration typically predicts steepening ra¬ 
dio spectra due to the decreasing ratio of acceleration 
efficiency to radiative losses for higher energy electrons. 
In order to confirm such models a detailed knowledge 
of the high frequency radio halo spectrum and morphol¬ 
ogy is desirable. The classical secondary injection mecha¬ 
nism is the decay of charged pions which were produced 
in hadronic interaction of relativistic protons with the 
background gas (Dennison 198C|, and others). However, 


Zeldovich 1980) of the hot electrons in the cluster plasma, 
as was also pointed out by Liang et al. (200C|). Our goal is 


to discuss the effects arising from the interference of these 
two radiation mechanisms. 

We use Hq = SOhsokms"^ Mpc“^, and assume Hq = 
65kms“^Mpc“^ in numerical examples. 

2. Thermal Sunyaev-Zeldovich effect 

The CMB has a perfect blackbody brightness spec¬ 
trum /cmb(j^) = ioi{x), with X = hi^/{kTcmh), *o = 
2 {kTcmhY/(h c)"^, and i{x) = x^l(eA — 1). The inverse 
Compton interaction of CMB photons with hot electrons 
in a galaxy cluster modifies the CMB spectrum. For small 
optical depth this is well described by 


= io g{x) y, with y = 


(7X 

TOe 


[ dl He kTe, 
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Fig. 1. Obs erved radio halo fluxes of the Coma cluster as 
compiled by Thierbach et al. ( 2002| ). The three radio halo 
emission models (A,B,C) described by Eqs. |[^ are dis¬ 
played in comparison to the isotropised synchrotron ker¬ 
nel (see text). The SZ decrement (a negative flux) is also 
shown. The solid angle over which the observed fluxes were 
integrated may vary from observation to observation. For 
the theoretical models, the flux from the full cluster vol¬ 
ume up to the assumed position of the accretion shock is 
given. 


the Compton y-parameter. rie, and Te are the electron 
density and temperature, respectively, I is the coordinate 
along the line-of-sight, and 

7^ - 4^ « -2 for x « 1, (2) 

describes the spectral distortionsjj 

The area integrated SZ distortion of the CMB by a 
galaxy cluster can be described as a SZ-flux (even if this 
is negative in the Rayleigh-Jeans part of the CMB): 


Fsz{i^) 


2 (Tt *0 Fth.e 
3 TOe 


9{x), 


(3) 


where Eth,e is the total energy of all thermal cluster elec¬ 
trons, and D is the angular cluster distance. 

In our examples, we assume the cluster to be isother¬ 
mal, spherically symmetric, and its electron density well 
described by a beta-model ne{r) = ne,o [1 + 
up to the accretion shock wave located at radius r = 
R. For Coma, we adopt kT^ = 8.2 keV, ne,o = 3 • 
10“^ cm“^ 


= Mpc =10', /3 = 0.8 (iBriel 


et al. 1 3 ^ ), R = 5 Mpc ( FnBlin et al. 1998| ) , and D = 

^ In the case of hot clusters (~ 10 keV) relativistic effects 
modify the shape of g{x). However, the corrections are neg¬ 
ligible on the Rayleigh-Jeans side of the CMB spec trum for 
kTe <C MeV (e.g. see Fig. 6 in Enfilin & Kaiser 2000), so that 
Eq. ^ is an accurate estimate in our frequency regime (x ^ 1). 
The relativistic electron population in Coma is also no t able 
to modify t his behaviour significantly ([Enfilin fc Kaiser 2000; 


Blasi et al, 2000; Shimon & Rephaeli 


2002D 



Fig. 2. Radio spectrum of the Coma cluster in v F{v) 
for better display (see Fig. |^. Solid (dotted) curves rep¬ 
resent the total cluster emission (CMB decrement) of the 
three halo emission models (A,B,C) combined with the SZ 
decrement. 


140 Mpc . This leads to = 7.7 • 10®^ erg, 

and a SZ-luminosity of 

Eir^{v) = 3.05(x) h-^iy « -1.9 • (4) 

where r'cHz = r^/GHz. The same parameters give a central 
y-parameter of y^odS = 0-85 • 10“"* which is in good 

agree ment with the measured value of = (0.7...1.0) ■ 


Herbig et al. 

1995; 

Mason et al. 

2001 

; De Petris et al. 


^ 0021 ) . 


3. Coma radio halo 


Cluster radio halos have steep radio spectra (Fhaio cx v~‘^) 
with spectral indexes s > 1. In the case of Coma, there 
seems to be a further steepening beyond 1.4 GHz. The 
current data as compiled by Thierbach et al. ( ^002 ) is dis¬ 
played in Fig. 1^ in comparison with three model spectra. 
A) A steep power law spectrum: 


= 0-723 ^GHz^Jy 

B) As A) but with a soft high frequency steepening: 

pComa/,,\ _ n .700 .,-1-34-0.lmax(0,lnj/GHa) t,. 

-^haio, 2 (i^l - b.723r'Qjj^ Jy 

G) A flat power law spectrum with exponential cutoff: 


(5) 


( 6 ) 


= 3-42 i^ghz® exp(-VJ^GHz) Jy- (7) 

The last spectrum is a fit to the data by Thierbaclj 


3t al. (2002), following a functional form proposed by 


Schlickeiser et al. (1987) as a model for a continously 
in-situ accelerated electron population with exponential 
high energy cutoff. It should be noted that the functional 
form of Eq. ^ was estimated with the monochromatic ap¬ 
proximation for the synchrotron emissivity of monoener- 
getic electrons. With the correct calculations, therefore. 
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Fig. 3. Radial brightness profile of our toy-model of the 
Coma cluster at different frequencies as it results from 
combining the radio halo emission and the SZ decrement 
profiles. The total emissivity follows the curve of model B 
displayed in Fig. ||. At 0.3 GHz (30 GHz) the practically 
unmodified radio halo (SZ decrement) profile is visible. 
The thin line gives the 5 GHz model for a = 1.1, [3 = 0.85, 
illustrating the strong dependence of the profiles on the 
parameters. 


the high frequency cutoff would be much smoother.^ In 
order to get a similar sharp cutoff in the emission spec¬ 
trum, a sharp cutoff in the electron spectrum is required, 
as e.g. in the time depen dent in-situ model investigated 
by Brunetti et al. (^001 ). Emission model G reproduces 
the observed data only if the SZ decrement is neglected 
(see Fig. ||). However, as explained above, this is not a 
problem for the continuous in-situ model, since correctly 
estimated, it produces a softer cutoff than given by model 
C. In model G the SZ-effect would dominate above 4.7 
GHz and lead to a decrement at the observed frequency 
of 4.9 GHz in contrast to the data. In model B (A) the 
change of sign of the total radio fluxes is at 5.7 GHz (6.2 
GHz). Model A predicts too much flux at high frequencies, 
but model B seems to be consistent with the data. Thus, 
only a weak spectral steepening of the radio halo emission 
above 1 GHz is required by the data. 

Such a weak steepening can also be explained in the 
secondary electron model. There the radio halo spectrum 
steepening may occur due to a steepening of the parent 
cosmic ray proton spectrum. 


4. Radial brightness profiles 

The various radio emission processes of galaxy clusters can 
partly be separated by their morphological difference. The 
SZ-emission has a spatially broad brightness distribution. 


since it is proportional (in an isothermal model) to the 
line-of-sight projected electron density. The radio emission 
seems to be as concentrated as the Bremsstrahlung X-ray 
emission of the thermal cluster gas: /halo oc where 
a « 1 (Govoni et al. 2001). Since the X-ray emissivity is 
proportional to rig, we assume the radio halo emissivity 
to exhibit the same scaling. Note, that secondary electron 
models would usually predict a steeper scaling (a « 1.3 
by Dolag & EnBlin (2000), but see Miniati et al. (2001)), 
which therefore seem to be disfavoured (Govoni et al. ^001 ; 
Brunetti 2(^2| ). 

In Fig. ^ we used a toy model to illustrate the conse¬ 
quences of the SZ decrement on the properties of the ra¬ 
dio brightness profiles as a function of frequency. For that 
purpose we assume a uniform emission spectrum given 
by model B (Eq. ^ and the above mentioned radial scal¬ 
ing.^. At low (^ GHz) frequencies the radio halo emission 
dominates. With increasing frequency a negative SZ-bowl 
surrounding the residual of the central visible radio halo 
emission appears and grows (dotted line in Fig. ^). At 10 
GHz the SZ decrement dominates at every angular posi¬ 
tion, but its central part still reveals the presence of the ra¬ 
dio halo emission by exhibiting a dip in the SZ decrement. 
At 30 GHz the radio halo emission is insignificant, so that 
this frequency can safely be used for SZ-measurements. 

Even at frequencies, at which the total radio emission 
is zero or negative (^ 5 GHz), a central positive emission 
region can exist if the frequency is not too high (< 10 
GHz). Thus, a positive flux may be attributed to such a 
radio map, if only the region with positive flux is included 
into the map analysis. Especially, since the negative bowl 
around the central region is relatively shallow (it domi¬ 
nates only the total flux due to its larger area) there is 
always the danger that it is regarded as noise, or as a 
remnant of imperfect map reconstruction in the case of an 
interferometric measurement without zero-spacing flux. In 
the case that the radio halo luminosity of a cluster is as¬ 
sumed to be equal to the flux within the positive region, 
the real halo luminosity would be underestimated (e.g. by 
34% at 5 GHz in our toy-model). 

In the frequency range of several GHz, there exist a 
relatively sharp outer edge of the positive brightness re¬ 
gion, which (in our toy-model) is located at 0 « 70' for 3 
GHz, and at 9 ze 30' for 5 GHz, whereas the synchrotron 
emission extends up to 6* = 120'. The exact location of the 
edge of the positive emission region depends sensitively on 
several of the model parameters (rc, /3, a, see Fig. ||) and 
on here not modelled details of the radio halo emission 
spectrum as a function of cluster radius. However, a sim¬ 
ilar frequency trend of the edge of positive emission can 
be seen by visual inspection of the Coma radio maps at 
1.4 GHz ( Deiss et al. 1997 ) which exhibits positive emis¬ 
sion up to 0 ^ 45', and the maps at 2.7 and 4.9 GHz 


^ The synchrotron spectrum of a monoenergetic electron 
population in an isotropic distribution of magnetic fields is 
displayed in Fig. A synchrotron spectrum can not be more 
strongly bent than this emission spectrum kernel. 


® The resulting spatially constant spectral index at low (< 
GHz) frequencies is not exactly what is observed since a ra- 
dial spectral steepening between 0.3 and 1.4 GHz is found by 


(Giovannini et al. 


1993). 
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( [Thierbach et al' ^002 ) which seem to be positive only 
up to 0 ~ 20' and 9 « 10', respectively. Although there 
is likely a contribution to this shrinking due to spectral 
steepening at larger radii, our estimate and the actual SZ 
measurements argue that a signihcant fraction of the emis¬ 
sion area shrinking is due to the unavoidable SZ-effect. 
The exact amount is difhcult to estimate due to our still 
limited knowledge of the exact radio halo emission pro¬ 
file. However, the prediction of the SZ-decrement profile 
is quiet robust for the central cluster regions due to the ob- 
servationally well constrained j/-parameter with an error 
of only ^ 10%. Therefore, it is possible to state that for 
the measurements of Thierbach et alj ( 2002| ) at 2.7 GHz 
and 4.9 GHz the central SZ-decrement is of the order of 
30% and 50% of the l-tr noise level (per beam as shown in 
their Fig. 2 and 3), respectively. Since the SZ decrement 
strength decreases radially only slowly, it should affect the 
outer regions of the rapidly decreasing radio halo^. Thus 
it should be statistically detectable there due to the suf¬ 
ficient large number of resolution elements covering the 
cluster outskirts. A careful analysis of high signal-to-noise 
high frequency maps of radio halos are expected to reveal 
an extended decrement at outer cluster regions. 


5. Conclusion 


Measurements of cluster radio halos at several GHz are 
contaminated by the thermal SZ effect, or vice versa, at 
~ 10 GHz, measurements of the SZ effect have to take the 
possible presence of cluster radio halos into account. The 
morphology of the radio emission of a galaxy cluster at 
such frequencies, where the total radio halo emission and 
the total SZ decrement roughly compensate each other, 
is that of a central emission region surrounded by a shal¬ 
low negative bowl. This complex morphology can lead to 
ambiguities in observational determinations of radio halo 
fluxes, since the result depends on the adopted procedure 
(flux from the full cluster area, or only flux from regions 
of positive brightness). We therefore stress the need to 
communicate up to which radius a cluster radio halo flux 
was integrated. It would be best, if such measurements 
would be provided for several radii. Radio halo emission 
measurements should carefully be corrected for the un¬ 
avoidable SZ decrement (see also Liang et al. 200C). 

Using the Coma cluster as an example we demon¬ 
strated that the existing radio maps at 2.7 and 4.9 GHz of 
the Coma radio halo should be affected by the SZ decre¬ 
ment. This is very likely part of the reason (in combination 
with radial spectral steepening of radio halo emissivity and 
a decreasing S/N ratio) why the apparent radio halo size 
decreases with increasing frequency. The observed sharp 
spectral steepening of the total radio flux of Coma seems 
to be partly - but not completely - due to the SZ effect. 


E.g. at the location of the S-cr brightness contour, located 


at 10' (7') in the 2.7 GHz (4.9 GHz) map of Thierbach et al, 
(2002), with a brightness of 0.038 (0.067) mjy/arcmin , the 


negative SZ contribution is —0.0026 (—0.011) mJy/arcmin^ 


This implies that only a weakly bended emission spectrum 
is required to explain the data, and consequently the con¬ 
strain arising from the apparent strong spectral steepen¬ 
ing on some proposed radio halo formation scenarios as 
continuous in-situ acceleration and secondary electron in¬ 
jection is relaxed. 
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